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In this paper we present a detailed and exact study of the probability density function P(F) of 
the total force F acting on a point particle belonging to a perturbed lattice of identical point sources 
of a power law pair interaction. The main results concern the large F tail of P(F) for which two 
cases are mainly distinguished: (i) Gaussian-like fast decreasing P(F) for lattice with perturbations 
forbidding any pair of particles to be found arbitrarily close to one each other; (ii) Levy-like power 
law decreasing P(F) when this possibility is instead permitted. It is important to note that in the 
second case the exponent of the power law tail of P(F) is the same for all perturbation (apart from 
very singular cases), and is in an one to one correspondence with the exponent characterizing the 
behavior of the pair interaction with the distance between the two particles. 



PACS numbers: 02.50.-r,05.40.-a,61.43.-j 



INTRODUCTION 



The knowledge of the statistical properties of the force 
acting on a particle belonging to a gas and exerted by 
all the other particles provides important information in 
many physical contexts and applications. Typical exam- 
ples are: (i) distribution of the gravitational force in a 
gas of masses in cosmological and stellar astrophysical 
applications Q, 0, 0j_ (ii) distribution of molecular and 
dipolar interactions Q in gas of particles, (iii) theory of 
defects in condensed matter physics 0, and (iv) granu- 
lar materials The first seminal work in this field 
was due to Chandrasekhar [lj and deals, among many 
subjects, with the gravitational force probability distri- 
bution in a homogeneous Poisson spatial distribution of 
identical particles. By studying the characteristic func- 
tion of the sum of the stochastic forces due to the single 
particles, the probability distribution of this total force 
is exactly found to be given by the so-called Holtzmark 
distribution, which is a three-dimensional analog of the 
one-dimensional fat tailed stable Levy distributions. In 
0,0, H,@ approximated extensions, to different branches 
of physics, of this approach can be found for more com- 
plex particle distribution (i.e., point processes) obtained 
by perturbing a homogeneous Poisson point process. In 
this paper we present a study of the total force probabil- 
ity distribution for a very different class of spatial par- 
ticle distributions (i.e., point processes), the perturbed 
lattices of point particles, in the case in which the pair 
interaction decays spatially as a general power law. We 
think that this study can be very useful for application in 
both solid state physics (e.g., in the case of Coulomb or 
dipolar pair interaction) [4( , and cosmology where n-body 
gravitational simulations (introduced to study the prob- 
lem of "structure formation" due to gravitational collapse 
from primordial cosmological mass density fluctuations) 
are performed usually starting from suitable perturbed 
lattice initial conditions 10]. We limit the study to the 
one-dimensional case in order to avoid difficulties related 
to the anisotropies of higher dimensional lattices. How- 



ever the exact results we present in Id are suggestive of 
the behavior of the same quantities in higher dimensions. 
In fact one can see [llj that the change of spatial dimen- 
sion only renders calculation not explicitly performable, 
keeping qualitatively the behavior we present below. 



II. DEFINITIONS AND FORMALISM 

In order to approach in the proper way the problem of 
the global force probability distribution in a perturbed 
lattice of particles interacting via a power law spatially 
decreasing pair interaction, let us consider firstly a gas of 
such identical particles with microscopic density 



n(x) = S(x — Xi) 



where Xi is the position of the i particle. Let us assume 
that the average number density uq = (n(x)) > (where 
the average (...) is to be intended as an ensemble average) 
is well defined (i.e., the particle distribution is uniform 
on sufficiently large scales 0). We then suppose that 
particles interact via a pair force f(x) depending on the 
mutual pair distance x as 



-c- 



|o+l 



This means that f(x) gives the force exerted by a par- 
ticle in the origin on another particle in x (the force is 
attractive if C > 0, and repulsive if C < 0). Therefore 
the force-field in the point x of the space will be 



T{x) 



la+l 



= C / dyn(y)- 



V 



\u- 



la+l 



(1) 

where the last integral is over all the space. However, 
from the last expression of Eq. JIJ, being no > 0, we have 
that for a given realization of the stochastic density field 
n(x), the infinite volume limit of T(x) is not univocally 
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defined for a < 1 (i.e. the integral in absolutely diverg- 
ing, and its value depend on how this limit is taken). The 
same feature is present in higher spatial dimensions. For 
example in d = 3 the same problem is present for a < 3 
and in general a < d in d dimensions. For instance this 
is the case of the gravitational force in a self-gravitating 
homogeneous gas of identical masses |jj , and of Coulomb 
interaction in the one component plasma (OCP) |l2( of 
identical electrical charges both in the disordered and 
the ordered (i.e., the Coulomb lattice 01) phases. This 
problem is well known in condensed matter physics about 
the OCP. However, in this case, the problem is automat- 
ically solved by the presence in the physical system of 
a uniform background charge density nb(x) = —no with 
opposite sign with respect to the identical charged parti- 
cles and such that to conserve global charge neutrality in 
the system. Once the attractive force of the background 
is considered on a charged particle together with the re- 
pulsive forces exerted by the other particles, the problem 
of the infinite volume limit of F is solved and its value 
unique. For what concerns the self-gravitating systems, 
in Newtonian gravitation an analog of the uniform back- 
ground of the OCP (i.e. a negative uniform mass density 
rib(x) — —ng such to generate a repulsive force on the 
particles) does not exist, and it has to be introduced in 
the system artificially to regularize t he p roblem (an ap- 
proach usually called Jeans' swindle 14]). However this 
negative background comes out naturally, as an effect of 
space expansion, when the gravitational motion of parti- 
cles is described, starting from the equations of general 
relativity, in comoving coordinates in a qua si- uniform ex- 
panding Einstein - De Sitter universe [15( which is the 
main model of universe used in cosmology. In practice 
considering the presence of such balancing background 
will give for T the following expression: 

F( X ) = cJdy5n(y) ]y y _-* +1 , (2) 

where Sn(x) = n(x) — no- This makes T to be defined 
also for smaller a depending on the small k behavior of 

the power spectrum S(k) ~ (\Sn(k)\ 2 ^ of the density 

field where 5n(k) is the Fourier transform of Sn(x). By 
studying the large distance scaling behavior of the inte- 
grated fluctuations of n(x) 0, it is simple to show that, 
assuming S(k) ~ k@ at small k, in d dimension T is a well 
defined statistical quantity (i.e., its value does not depend 
on the way in which the infinite volume limit is taken) 
for a > (d - 0)/2 if j3 < 1, and a > (d - l)/2 if (3 > 1 
(see also the discussion in Appendix II on the definite- 
ness of the force T respectively in the shuffled lattice and 
the homogeneous Poisson particle distributions). Note fi- 
nally that taking the infinite volume limit symmetrically 
with respect to the point x on which the force is calcu- 
lated the background gives a zero net force on the point 
x. Therefore the value of T obtained calculating Eq. (JTJ 
taking the infinite volume limit symmetrically with re- 
spect to the point x gives automatically the well defined 



value obtained by Eq. (i.e., subtracting the effect of 
the background). This preliminary discussion of the sta- 
tistical definiteness of the force T is useful to justify the 
symmetrical way in which we take the infinite volume 
limit in the one-dimensional shuffled lattice case we an- 
alyze in the rest of the paper. The statistical properties 
of T we will find in this peculiar way (no background 
and symmetrical limit) coincide with those of the case in 
which the effect of a negative background is considered 
independently of the way in which the infinite volume 
limit is taken, which therefore can be considered the real 
physical case. Moreover, from the above considerations 
we deduce that, as in the shuffled lattice S(k) ~ k 2 at 
small k (see ^fj)> the results are valid for all values a > 
in d = 1. 

Let us take, therefore, a Id regular chain of 2N + 1 
unitary mass particles with a lattice spacing a > (we 
will take eventually the limit N — > +oo), i.e., the position 
of the n th particle is X n = na. Therefore the microscopic 
density can be written as 

N 

n in {x) = 2_j - na ) ■ 
n=-N 

Clearly the average density of particles in the system 
is no — 1/a. We now apply an uncorrelated displace- 
ment field (i.e., a random shuffling) to this system, i.e., 
a random displacement U n is applied to the generic n th 
particle independently of the other particles. This dis- 
placement field is completely characterized by the one- 
displacement probability density function (PDF) p(u) 
(i.e., Prob(u < U n < u + du) = p{u)du). After the appli- 
cation of the displacements the new microscopic density 
will be: 

N 

n(x) = 5(x — na — u n ) , (3) 

n=-N 

the u n 's being the realizations of the random variables 
U n all extracted from p(u) independently one of each 
other. We consider the case p(u) — p{—u) for simplicity. 
Equation (0 says that the particle originally in X n = na, 
after the displacement will be in X n = na + U n . For the 
analysis of spatial density correlations in such a system 
see [l|| . Let us call q n (x) the PDF of the position of the 
n th particle (i.e., Prob(x < x n < x + dx) = q n (x)dx). 
Clearly it is given by 

q n (x) = p(x - na) . 

Let us now assume, as above, that the n th particle 
creates a force-field in the point x of the type: 



with a > and C a constant. Therefore the total 
stochastic field J- generated at a generic point x of the 
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space by all the system particles is: 

N 



T(x) = C 



n=-N 



X n - X 
\X n - 



(4) 



Note that it is a sum of random variables. Let us call 
W (F) its PDF; it will be given by 



W (F) 




It is immediate to see that T{x) is the sum of independent 
random variables C(X n — x)/\X n — x\ a+1 . However, as 
the PDF's q n (x) change with n, these variable are not 
identically distributed. As shown below, this, together 
with the fact that T needs not a normalization in N 
to be well defined in the large N limit, are the reasons 
why we do not obtain in general an exact Gaussian or 
Levy limit 0,0] for the Wq{F). In order to study the 
asymptotic behavior in F and N, it is usual to introduce 
the so-called characteristic function of T , i.e., the Fourier 
transform (FT) of W (F): 



W (k)= I dFW (F)e lkF 

N 



f°° ( y — x 

1^ dyp{y ~ na) 6XP \ iCk \y-x\ a + 1 



(5) 



By studying the small k behavior of the single integrals 
in Eq. and taking appropriately the limit N — > oo 
we can deduce the moments and the large F behavior of 
W (F). 

However, we will study a slightly different and more 
difficult problem, which is of particular interest if we want 
to study the dynamics of the system particles under the 
effect of only this mutual force. We study directly the 
statistical properties of the stochastic force acting on one 
generic system particle. In particular we calculate the 
total force acting on the particle initially located at the 
origin of the space and displaced to Xq = Uq: 

-N,N —N,N 



n#0 



\X n -X \ a+1 



In this way, taking the limit N — > oo, we get the sym- 
metric infinite volume limit with no negative background 
of Eq.lpQ) with x — Xq, which is, as explained in the first 
paragraph of this section, equal to Eq.J2J with arbitrary 
way of taking the infinite volume limit, and therefore giv- 
ing the force on any particle belonging to the system with 
the uniform balancing background. Note that now, be- 
cause of the presence of the variable Xq in each term of 
the sum @, J 7 is no more a sum of independent terms. 



However we show below how to reduce the problem to 
that of a sum of independent stochastic terms, by in- 
troducing the concept of conditional probability density 
function. The solution of this conditional problem will 
give also the way to face the study of the first uncondi- 
tional case given by Eq. Q . 

About T in Eq. J§}, we want again to find the PDF 
W(F) of the value F of this force. As before, since the 
displacements applied to the particles are independent 
one of each other, we have the exact relation: 



+ OC 




A' 



Y[ dx n p(x n - na) 



N,N 



-N 



X n Xq 



n/0 



\x n - Xq\ 



a+1 



that, through a simple change of variables A„ = x n — xq 
or n 7^ 0, can be rewritten as 



/+oo 
dxQ P (xo) 



— N.N 

Y\ dA n p(A n + xq - na) 



—N,N 



Let us analyze the behavior of the conditional PDF 
P(F;xo), conditioned to the fact that the particle on 
which the force is evaluated is at Xq = xq\ 



P(F;x ) 



-N,N 

Yl dA n p(A n +x - na) 



-N.N 



;S\F-C £ 



A, 



(7) 



n^O 



In this way, once xq is fixed, the total force J- is the sum 
of independent contributions 



fn = C- 



(8) 



It is interesting to see in which cases T satisfies the cen- 
tral limit theorem. We will see that it never satisfies this 
theorem even when its PDF is rapidly decreasing at large 
values. More precisely, we will see that even in this last 
case its PDF is dependent on the details of p(u). 

For the sake of simplicity of notation let us assume 
in the rest of the paper that C = 1 (the repulsive case 
C = — 1 will be trivially deduced from this). By per- 
forming the simple change of variables given by Eq. JSJ , 
it is possible to find the conditional (i.e., conditioned to 
^0 = ^o) PDF g n (f;xo) of the single stochastic force 
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generated by the particle in X n on the particle fixed in 

9n(f; xo) = P ( + x ° ~ na ) • ( 9 ) 

Clearly also the forces /„ so distributed are independent 
of one each other. The support of g n {f', Xq) can be simply 
deduced from the one of p(u). 



III. LARGE F ANALYSIS AND LIMIT 
THEOREMS 

Let us take the FT of Eq. Q in order to evaluate 
the conditional characteristic function of the stochastic 
variable F: 



P(k;x ) 

—N,N 

- n 

n/0 



+00 



dFP{F;x )e 



IkF 



(10) 



dx p(x + xn — na) exp ik-, — : — -r 
1 \x\ a+1 



Note that the quantity 



dxp(x + xq — na) exp ik 



(11) 



df g n {f;x )e 



ik f = f e ik f 



\ e ;n;i|) 



where g n (f', xo), given by Eq. JJjJ, is the conditional PDF 
of the field / felt by the particle in xq due to the particle 
in x n , is the conditional characteristic function of this 
field /. Moreover {a{f)} n - Xo = f_™ df g n (f;x )a(f). It 
is important to note that, if x$ = (i.e., the particle 
on which we calculate the force is stuck at the origin), 
the condition p(u) = p(—u) on the displacements of any 
particle would imply 

+00 



dx p(x + na) exp ik- — ; — ; 

» V M" 4 

f + °° ( x 

= / dxpix — na) exp — ik-. — ; — rr 

Consequently, if one fixes x = Eq. (fTT7|> becomes: 
W(k) = Yl J dxp(x - na) exp |^fc ^ a+1 J ■ 

However the shift xo ^ of the particle initially at the 
origin, and on which we calculate the force, breaks this 
symmetry, which is anyway recovered when the aver- 
age over p{xq) is performed to proceed from P(F; xq) to 
W(F) (however we will see that this is a further source 
of noise when we calculate explicitly the variance of the 
force F). 

In order to proceed into the analysis of the PDF of T, 
we have to distinguish two basically different cases: 



1. Non-Overlapping condition (NOC): No parti- 
cle can be found arbitrarily close to any other par- 
ticle; i.e., the supports |20l| respectively of p(u) and 
of p(u — na), for all integer n 0, have an empty 
overlap. The main case of physical interest in this 
class of displacement fields is when 30 < uq < a/2 
such that p(u) — for \u\ > uq', 

2. Overlapping condition (OC): Particles can cross 
one each other and at least one pair of particles can 
be found arbitrarily close to one each other; i.e., the 
supports respectively of p(u) and of p(u — na), for 
at least an integer n ^ 0, have a non-zero overlap. 
The main case of physical interest in which this 
happens is when Be > such that p(u) > for all 
M < a/2 + e. 

We will see that in the first case we obtain a rapidly de- 
creasing W(F) even though there is no constraint toward 
Gaussianity in the large N limit, while in the second case 
we have a power law tailed W(F) similarly to that of the 
three-dimensional Holtzmark distribution [J. 



IV. DETAILED ANALYSIS OF EQ. (fTTTl 

Let us analyze the single factor of Eq. (|10|) which, 
as aforementioned, is the FT of the conditional PDF 
g-n(f',xo) of the force felt by the particle in Xq = xq 
due to only the particle in X n : 

\ dx nix + xq — na) exp ik- — : — r-r I 
J -00 V \x\ a+1 J 



= ( exp I ik 



(12) 



where (s(x)) denotes the average of the function 
s(x) over the shifted PDF h nxo (x) = p{x + a; — 
na). In practice, if we indicate with simply (s(u)) — 
dup(u)s(u), then we can say that 

( s ( x ))n, Xo = (s(u + na-x )) ■ 

We want to study Eq. Ijllfl in the limit of small k. Sim- 
ilarly to what pointed out in the previous section, the 

small k behavior of (exp ^fc ^S+i ^y is different in 

the two cases in which, as a consequence of displace- 
ments, the pair of particles initially in x = and x = na 
cannot or can be found arbitrarily close to one each other, 
i.e., respectively if the supports of p(u) and p(u—na) have 
an empty or a non-zero overlap. 

Let us start with the case (i). If 3 < u* < \n\a/2 such 

that p(u) — for |tt| > u* , the exponent of exp (ik \ x \L+i \ 

can take only limited values in the integral l|12|) (i.e., the 
support of g n {f',Xo) is restricted to only a finite interval 
of values of /). Note that in the given hypothesis (i), if 
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n > the quantity x can take only strictly positive val- 
ues, while if n < it takes only strictly negative values. 
In this case, if n > 0, we can write: 



exp ik 



\a+l 



— oo 



, . , (ikx- a y 

dxpix + xq — na) > ; — 

' ml 



m—Q 



m=0 



z — / m \ 

where 

((»i + Tia-io) _am )= / dup{u){u-x + na)- am . 
J —u* 

If instead n < 0, Eq. l |13J t becomes 



+00 



{-iky 



n,x m=0 



m I 



((-x)- am ) 



+00 

E 

m=0 



(-iky 



■({-u-na + x o y am ) , 



(14) 



Note that, as p(u) = p(—u), we have 
{(-u-na + x )- am ) = {(u-na + x )- am ). In 
both cases 



1 0+1 



< (no + 



2m*)" 



for any m > 0, and therefore the series in Eq. ab- 
solutely converges. It is very important to note that, if 
u* < a/2 (i.e., no pair of particles can be found arbitrar- 
ily close to one each other) all the factors in Eq. (|10fl can 
be represented as a Taylor series 1)130 to all orders m. 
As shown below in more detail, this implies that when 
u* < a/2, W(F) has all finite moments and therefore is 
rapidly decreasing at large F. 

In the second case (ii) instead 30 < e < a/2 such that, 
at least Vit satisfying \n\a/2 — e < \u\ < \n\a/2 + e, we 
have p(u) > 0. In this case the particle initially at ±na 
and the particle initially at the origin can be found ar- 
bitrarily close. This implies that the quantity x (i.e., 
u — xq + na) in (|12|) is permitted to take arbitrary small 
values up to zero, and therefore in the last expression 
of Eq. Ill.'Sfl there would be an infinite number of diverg- 
ing terms of the last series. In other words the Taylor 
series sum in the second expression of Eq. <|13|) cannot 
be exchanged with the average operation (...}„ Xg , and 
we expect a singular part in the small k expansion of 

the average (^exp {ik ^prr J j ■ In order to find it, 



WO 



rewrite it as in Eq. (|11|) : 



6XP V k W+ 1 ) / ^ J d f 9n{f]x )e ikf 

+ °° fx 
dx p(x + xq — na) exp ik- — ; — — 

-oo V \*\ a+1 



Note that in this case, differently from the previous one, 
the support of <?«(/; £o) includes arbitrarily large values 
of |/| for which, using Eq. @, we have 



= P(x -na) m _ {a+1)/a + o(|/r (a +1)/a) . 
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Let us call M = [a x ] the integer part of a 1 . By using 
the results presented in Appendix I we can conclude that 



exp ik-r—. — -r 



(15) 



n;XQ 



M 



= E ^T- (( u + na - x o)~ am ) + S n (k; x ) , 
* — ' ml 

m— 

where S n (k;xo) contains all the terms of order higher 
than M, including the singular part of the small k ex- 
pansion of (exp (ik j^ja+i ^ ^> which is of order 1/a in 
k. By using Eq. (|3"2*|) . we can finally write: 



S„(fc;^ ) 



(16) 



T (xo r"! } , fc 1 /" + oik 1 '*) for odd M 

ar[(a+l)/a]cos( ° f M ir) 

(-l) M/2 «P(.Jo-_n a ) k l/ a + o(fc l/o) for cvcn M 

ar[(a+l)/a]sin( ° 2 V ) 

Here, for simplicity, we have excluded the case in which 
exactly M = 1/a for which we have logarithmic correc- 
tions in k to the above equations. 



V. FINDING W(F) 

At this point we can go further and classify the possible 
behaviors of P(F; xq) and W(F). 

Basically we again distinguish the following two cases: 

1. 3e > such that V|u| > a/2 — e one has p(u) = 0; 

2. 30 < e < a/2 such that, at least Vu satisfying a/2 — 
e < M < a/2 + e, p(u) > 0. 



A. Case 1: fast decreasing W(F) 

In this case the system satisfies the NOC and all the 
factors in Eq. i|10|) can be expanded in the Taylor series 
(O and {nj for all different n. 
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We can then write 

r + co 

P(k;x) = / dFP(F;x)e tkF 



(17) 



JV 



] J (c 



n 

n=l 
+00 

E 



E 

,m=0 
Z! 



(z/fc) r 



((ji-i + m)-"™} 



({u + x + na)- al ) 



where in the average (...) over the displacement u we 
have used the symmetry property p(u) = p(—u). Note 
that from Eq. (|17|) . we have P(k; —x) = P'(k;x), where 
indicates the complex conjugate of A. 
By calling again u* < a/2 the maximal permitted dis- 
placement for each particle (i.e., the support of p(u) is 
included in [-u*,u*]), we can find W(k) = T[W(F)] by 
simply calculating the following average 



W(k) = / dxp(x)P(k;x) 



(18) 



It is simple to verify that, as P(k;—x) = P^(k;x) and 
p(u) — p(—u), the function W(k) is real and W(k) = 
W(— k). The Taylor expansion in k of W(k) is obtainable 
from Eas. i|17|) and |JTSJ|. Since it is a real function and 
is the characteristic function only even powers of k are 
present. In particular the coefficient of the k 2 term is 
—T 2 /2 where 



h(T) 



+ OG 



dF W(F)h(F) . 



Actually, rigorously speaking, we should show that all 
the coefficients of the Taylor expansion of W(k) are con- 
vergent to finite values in the limit N —> 00. It is simple 
to show it by expanding the terms ((u ± xq + na)~ am ) 
of Eq. (|17|) in Taylor series of (it ± xo)/na for n > 1 
which is justified by the fact that in the given hypothesis 
\u\ + \xq\ < 2u* < a, and considering that 



+00 

E 

m=0 



(ik) 7 



■(na)~ 



^ (-ik) 1 

^ l\ 

1=0 



(na)~ 



= 1 , Vn > 1 . 



Therefore we conclude that we can write W(k) in the 
following form: 



W(k) 



+ OO 



^ n (2ny. k2n 

n=0 y ' 



It is simple to see that 

JV 



^ = E {«(«- »+~)- Sto U B as) 

n=l 

-{((u-x + na)- a ) u {(u + x + nar«) u ) x } 



l,N 
n<l 



x + na)- a ) u -((u + x + na)- a )J 

x [((u-x + la)- a ) u -((u + x + la)- a ) u ]) x , 
where for clarity we have redefined 

' ( a ( u )) u = I-u' dua(u)p(u) 
' ( a ( x )) x = f- u * dxa(x)p(x) 
. ((K u ^)) u ) x du dx b(u, x)p(u)p(x) . 

It is matter of simple algebra to show that, for p(u) 
p(—u), Eq. (|T§)) can be rewritten as 



T 2 



N 



ir=E«( u - x + na ^ 2a )u - d u - x + na r a )u 



n=l 



1,N 



+ (((« ~x + na)- a ) u [((« - x + la)- a ) u 

It. I 

-((u + x + lar a ) u }) x (20) 

We see that the force variance is composed of two differ- 
ent contributions: the former, given by the first sum in 
Eq. i(2U)l , is mainly due to the fluctuations in the displace- 
ments u of all the sources of the force (in this term the 
average over x is only a smoothing operation), while the 
latter, given by the second sum, is determined basically 
by the fluctuations created by the stochastic displace- 
ment x of the particle initially in the origin on which we 
evaluate the force (in this term is the averages over u to 
play a role of simple smoothing) . 

It is interesting and useful in applications to calculate 
all the above expressions by evaluating all the terms in 
the sums in Eq. ij2U|) to the second order in (u ± x)/na. 
In order to do this, we use the following second order 
Taylor expansion for B <C A: 



B 

1+ A 



(A + B)—i =A-~< 

B 7(7+1) 



A -7 



1-7^ + 



From this, substituting respectively A with na and B 
with u ± x, we have that 



(u ± x + na\ 



(na) 



'7- 



u±x 7(7+1) /llil 



7 



Moreover we have that (((u ± x ) 2n+1 ) u ) x — f° r an Y 
integer n due to the symmetry p(u) = p(—u), while we 
have that (((u ± a;) 2 )^) = 2er 2 where a 2 = (u 2 ) u = 

f_ . duu 2 p(u) is the variance of the single displacement. 
Therefore we can write 



(na)" 



«( U ±x + na)^)J 
and 

(((u±x + na)- a )l) 

Henceforth 

({(u±x + na)- 2a ) u -(( 



1+7(7+1) 



(na) 2 



'. (na)' 



2d 



1 +a(3a + 2) 



(na) 



u ± x ■ 



na)-°)l) 3 



( na )2( Q +l) ' 

Moreover 

(((n-x + na)- a > u ((u± a; + Zar a >J :c 
~ (na)~ a {la)- a 

11 



1 + 



+ a(a + llcr" 



(la)(na) \(na) 2 (la) 2 

from which 

{{(u-x + na)- a )J{(u-x + lar a ) u 

~((u + x + la)- a )J) 



(21) 



2a 2 a 



2^-2 



(la)«+ 1 H Q + 1 ■ 

Using all these results in all the terms of Eq. i|2U[) . we 
obtain 

-1 2 " 



2 



2 2 
a a 



+ 2 



JV 

E 



1 



■i (na)2(«+ 1 ) ^ (na) a+1 

= 1 v y Ln=l V 7 

(22) 

It is simple to verify that both sums in Eq. (|22|l are con- 
verging for N — ► +00 for all a > 0, for which we can then 
rewrite 



El 



a 2 a 2 

2 2(a+l) 



[C(2a + 2) + 2( 2 (a + l) 



(23) 



where is the Riemann zeta function (note that for 
t — > 1 + wee have ((t) ~ l/(t — 1)). Again the first term 
is due to the fluctuations in the position of the sources, 
while the second one is due to the fluctuations in the 
position of the particle on which we are calculating the 
force. In particular In Eq. 12211 the generic term of the 
first sum give the relative weight of the n th nearest neigh- 
bor particles in determining the force on the particle in 
Xq. At last we can say that, in the case of displace- 
ments limited within a box well contained in a unitary 
cell around the initial lattice position, we can approxi- 
mate W(F) with a Gaussian PDF with zero mean and 
variance given by Eq. I|23l) . However, as already pointed 
out, there is no constraint, in the limit N — * 00, toward 
rigorous Gaussianity and non-Gaussian corrections are in 
general present. 



B. Case 2: power law tailed W(F) 

As shown above, this is the case in which the OC is sat- 
isfied, i.e., particles are permitted to jump out of their 
initial lattice positions beyond the limit of the related 
unitary cell in such a way to be found arbitrarily close 
to some other particle. Note that this is always the case 
when the support of p(u) is unlimited, i.e., if p(u) > 0, 
Vu £ M. However the same kind of W(F) is also ob- 
tained if 3u* > a/2 such that p(u) > 0, Vu € [— u*,u*} 
and zero outside. The difference between these two sub- 
cases is only in the amplitude of the power law tail of 
W(F) but not in its exponent. In general, if the particle 
initially at the lattice site x = na is permitted, through 
displacements, to be found arbitrarily close to the par- 
ticle initially at x = 0, it will contribute to the product 
(|ll)fl through a factor of the type H15|) . If instead this is 
not permitted, it will contribute to (|10f> through a fac- 
tor of the form l|13(l or (|14J) depending respectively on 
whether n > or n < 0. In any case if M — [1/a], in 
order to find the main terms of the small k expansion of 
W(k) (so to determine the large F tail of W(F)), it is 
sufficient to truncate all the small k expansion of the dif- 
ferent factors in Eq. H10JI at most to the order M+l. For 
the sake of simplicity, let us limit the discussion to the 
case in which strictly aM < 1 in such a way to exclude 
logarithmic corrections in k. We can write 



P(k;x ) 

(OC) 
n 

(noc) 



M 

E 



(ik) 1 



n 



^(ikY 



(24) 

A(n, x 0l a)k 1/a 



\a+l 



l.Xn 



where A(n, xq, a) is the coefficient of the term ~ A; 1 /" in 
Eq. (|16l) , and the first product on n is on the particles in 
X n with n =/= which can be found arbitrarily near to the 
particle in X (i.e., satisfying the OC with respect to the 
particle in Xq), while the product on I is on the particles 
in Xi with I 7^ which have a positive minimal distance 
to the same particle (i.e., satisfying the NOC with respect 
to the particle in Xo). If p(u) > Vu S M all the system 
particles with n ^ are included in the first product. If 
instead p(u) > for u 6 [— u*,u*] with u* > a/2 and 
zero outside, the first product include only contributions 
from the particles with — 2u*/a < n < 2u* /a and n^0, 
while the others are included in the second product. The 
large F behavior of P(F; xo) and consequently of W(F) 
is completely determined by the singular term of order 
k 1 / in the small k expansion of P(k; xq). It is simple to 
see that up to the order k x l a 



P(k; x ) = c m (x )k m + c 1/a (x )k 1/a , 

where the the coefficients c m (xo) can be deduced 
by counting from Eq. (|24fl (in particular c m (xo) = 



8 



i m J rm (x )/m\, where F m {x ) = dF P(F; x )F m is 
in "' moment of P(F; xq) and m < M) and 



the -~> th 



Cl/c 



E 

-2u*/a<n<2u*/a 



A(n;x ) . 



where the formula includes also the case u* — > oo. The 
small fc expansion up to the order fc 1 /" of W(F) will be 
consequently 



VF(fc) 



M 

E & « 



where 
b m = 

bi/ a = 



dx p(x )c m (x ) 



-T m with m<M 



dxop(x )c 1/a (xo) , 



where W^Fji^' 

rf^oP(a ; o)-F m (a;o) = (^J 7m (x )'^. It is possible 

to evaluate explicitly by using Eq. ifTfijl . 

We are now in the situation to connect the singular 
term b\j a k x l a of the small k expansion of W(k) to the 
large F tail of W(F) by using directly the arguments in 
Appendix I. This gives simply 



W(F) ~ BF- l - 1/a 



with 



B = - 

a 



+00 



dx p(x ) 



E p ( x ° 

-2u* /a<n<2ii* / a 



•na). (25) 



Note that if the support of p(xo — na) is much larger 
than a and p(u) is smooth (i.e., approximately constant) 
on the scale a, we can approximate Eq. I|25[) with 



B = 



1 -p{0)a 



na 



Note that this last approximated expression is not de- 
pendent on the details of p(u) for u ^ 0. Finally, we 
can observe that we have obtained a power law tailed 
W(F) characterized by the same exponent of the case of 
a homogeneous Poisson particle distribution presented in 
Appendix II. The only differences are the two following: 

• The amplitude of this power law tail is reduced 
in the shuffled lattice with respect to that of the 
Poisson particle distribution, given by Eq. IpTTj) . of 
a factor 



+00 



n#0 

dxop(xo) ap(xQ — na) 

— 2u* / a<n<2u* / a 



l-p(0) a. 



• In the shuffled lattice we have this power law tail for 
each n > 0, while in the Poisson case the problem 
is not well defined for n < 1/2 (see Appendix II). 



VI. CONCLUSIONS 

We have presented a detailed study of the PDF W(F) 
of the stochastic force T generated by a randomly per- 
turbed lattice of sources of a scale invariant attractive 
pair interaction field f(x) = — Cx/\x\ a+1 with n > at 
distance x from the source. 

In general we distinguish two cases: 

1. The NOC is satisfied and no pair of particles can 
be found at an arbitrarily small reciprocal distance; 

2. The OC is satisfied and it exists at least one of such 
pairs of particles. 

In the first case we have a fast decreasing W(F) similar 
to a Gaussian PDF at large F, even though no constraint 
toward an exact Gaussian central limit theorem is found. 
In the second case a power law tailed W(F) is found. 
The unique exponent of such power law is directly related 
to the pair interaction exponent n, while its amplitude 
depends also on the lattice spacing a (with respect to 
the unit distance through which we measure x in f(x)) 
and in general on the shape of the perturbations PDF 
p(u). In particular in this case W(F) has a power law tail 
with the same exponent as the stable Levy distribution 
found in the Poisson case (see Appendix II) but with a 
reduced amplitude, even though, analogously to the case 
(i), no constraint has been found toward the stable Levy 
distribution. 

Some further general considerations have can now be 
done: 

• In the case in which the probability of finding arbi- 
trarily close to one each other, the large F behavior 
of W(F) is basically determined by the small x be- 
havior of fix) and not at all by the the fact if it 
is long range or not. Therefore if we considered a 
fast decreasing f(x) but with the same divergence 
in x = we would have deduced the same con- 
clusions about the exponent of the large F tail of 
W(F) [13. 

• For this reason, even if we consider a lattice per- 
turbed by correlated displacements, we expect to 
obtain the partition into the two cases (i) and (ii) 
above considered depending on the possibility or 
not to find pair of particles arbitrarily close. 

• Actually different cases for the large F tail of W(F) 
between the Gaussian-like "fast decreasing" and 
Levy-like "power law" tailed PDF with exponent 
(3 = (n + l)/n are possible in very particular cases. 
These cases correspond to the choice of p(u) such 
that p(u) > exactly for u £ [— a/2, a/2] and zero 
outside. By changing the limit behaviors of p{u) 
when u — > ±a/2 we can obtain different large F 
behaviors of W(F). In particular if p(a/2) > and 
finite (we consider p(u) — p(—u)) we have the same 
case as (ii) described above with (3 — (n+ I) /a. If 
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instead p(a/2) = 0, depending on the behavior of 
p(u) for u — > (a/2) - we will have different values of 
(3 but in general larger than (a + l)/a. If, finally, 
p(a/2) = +oo (in such a way that p(u) remains any- 
way integrable) in general we obtain [3 smaller than 
(a + 1)/ a (but always > 1 so that W(F) remains 
integrable). 



Appendix I: Fourier transform of power law tailed 
PDFs 



We are interested in the small k behavior of the char- 
acteristic function f(k) of a given power law tailed PDF 
f(x) which for large |a;| behaves as with a > 1. 

Let us call [a] = n > 1 the integer part of a. In this 
hypothesis f(k) has a regular Taylor expansion up to the 
order n — 1 followed by a singular term proportional to 
k a ~ x : 



/+oo n—l 
dxf(x)e ikx = J2 ■ 
-°° m=0 



(ik) m 



x m + f s (k), (26) 



where x m — f_^dxx m p(x) and f s (k) contains the sin- 
gular part of f(k) and at small k is an infinitesimal of 
order a — 1 in k (if a is an integer it contains also loga- 
rithmic corrections). 

Now we show that effectively at sufficiently small k, 
fs(k) ~ B/c Q-1 (where now a(k) ~ b(k) means that 
linife—,0 [^(fc) /h(k)] = 1) giving an explicit expression for 
B as a function of both the amplitude A and the expo- 
nent a. 

First of all, let us study the case of a function h(x) 
that can be written as 



h(x) = B\x\~P +h (x) 



(27) 



where B > 0, < (3 < I and ho(x) is a smooth func- 
tion, integrable in x = and such that x^ho(x) — ► 
for |a;| — > oo. This means that h{x) presents an even 
power law tail. In this case the small k behavior of 
h(k) = f_ °° dx h(x)e lkx is completely determined by the 



Fourier transform of B\ 



i.e. 



h(k) ~B dx\x\- e ikx 



In order to perform this Fourier transform we introduce 
the integral representation: 



-0 



1 



r(/3) Jo 



0-l-\x\z 



dzz p - L e 



(28) 



where T(f3) is the Euler Gamma function. Using Eq. I|28|) 
we can write: 

dzz 13 ' 1 I dxe lkx - yAz 



/-\-oo 1 POO p 



-oo 

2fc' 3-1 



dq-_ 



r(/3) Jo ^l + q- 



Using the general relation valid for < j3 < 1 



X dqT ^ = (29) 



! + ? 2 2cos(f^)' 



we can conclude 



h{k) = 



Bn 



T(/3) cos (§tt) 



fc/3-i + oik ' 1 ) . (30) 



If instead of using the integral representation l|28(l one 
used 



-P 



1 



one should obtain the alternative expression containing 
only Gamma functions: 



h{k) 



2"r(§) 



(^ -1 ). 



Note that the coefficient of the term is real and 

positive. Another important case is when the function 
h[x) has an odd non-integrable power law tail, i.e.: 



h{x) = B\x\-P[20{x) - 1] + h (x) . 



(31) 



where 9(x) is the usual Heaviside step function, B > 0, 
< j3 < 1, and ho(x) the same of Eq. (|27|l . By using the 
same integral transformation leading to Eq. (|30|l . we in 
this case we obtain: 



h{k) 



Bit 



T{f3) sin(|7r 



Ak^ 1 ). 



At this point we can go back to the problem of finding 
the dominant small k contribution of the term f s (k) in 
Eq. (|2"6*|l for the PDF f(x) decaying at large |x| as A|a;| -Q! . 
Note that now we cannot apply directly the argument we 
have used for the above function h(x). In fact if, from 
one side, also in this case we can write 



f(x)=A\x\- 



fo(x) 



with \x\ a fo(x) — » for \x\ — > oo, from the other side a > 
1 (for definiteness of probability) and fo(x) contains a 
non integrable singularity at x — such that to cancel the 
non integrable contribution of the A|cc| -a term at small 
x. In order to circumvent this difficulty we introduce the 
function 

g(x) = x n f(x) , 

where n is the integer part of a. In this way g(x) is 
similar to the function h(x) of Eq. (|27|l if n is even and 
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to the h(x) of Eq. (|31[1 if n is odd. Therefore, by defining 



as usual g(k) = /_ °° dx g( 



x)e 



ikx 



we can say that 



m 



An 



T(a - n) cos (^tt) 
if n is even, and 

An 



<a — n — 1 



,a — n — 1 



o(fc° 



r(a-n)sin(^»7r)' 



o(k 



a— n— 1 



if n is odd. Now in order to find the singular part f s {k) 
of f(k) it is sufficient to integrate n times g(k) (the inte- 
gration constants giving rise to the finite moments terms 
of f(k) in Eq. (|26|l ). In this way we obtain 



fs(k) 



(32) 



^ 7 \ k a 1 + 1 ) for even n 

r(a) cos! 7r ) v 7 



(-1)^+D/»A tt fca -l 
r(a) sin( a 2 n 7tJ 



o(fc 



a-l\ 



for odd n . 



where we have used the following property of Gamma 
function: r(x + l) = xT(x), hence [(a — l)...(a — n)T(a — 
n)] = T(a). Note that in both case the coefficient of the 
term fc"^ 1 is real. 



Appendix II: Force PDF in a Poisson particle 
distribution 



Let us consider the case in which the particles are dis- 
tributed on the line interval (—L/2, L/2] of length L fol- 
lowing a spatially stationary Poisson process with aver- 
age density p > 0. We want to know the PDF Wp(F) 
of the field 



F 



N 

E — 

i=i 



X; 



a+1 



(33) 



generated at the origin of the space by all the N system 
particles (we can consider N = pqL as the fluctuations of 
order \fp^L 1 due to the Poisson statistics, are completely 
unimportant for this problem in the large L limit). We 
will follow the procedure to find Wp(F) in three dimen- 
sions used by Chandrasekhar in l] for the gravitational 
force. Note that, as the positions of different particles are 
uncorrelated, the joint PDF j)(a;i, ...,xn) of the positions 
of the N system particles is simply: 



N 



p N (xi, ...,xn) 



rip( a 

i=l 



where p(xi) = l/L. Therefore we can write 



W P (F) = 



L/2 
-L/2 



' N 

n 



dxi 



N 



(34) 



Let us now study the characteristic function 

/+oo 
dFW P {F)e lkF . 
-oo 

By taking the FT of Eq. (O we obtain 

" ' L/2 dx 



Wp(k) 



L/2 



, ■ exp ik-r—. -r 

L \ \x\ a+1 



N 



By adding and subtracting 1 inside the square brackets, 
and taking the limit L — > oo with iV = p^L we arrive at 
the final expression: 



Wp(k) = exp 
Note that 



k x ' a 



-Po 



dt I 1 — exp i 



+ oo 



dt I 1 — exp i 



t 



\t\ a+1 



dt (1 - cos(^ Q )) 



tftttt~ 1_1 / a (l - cosu) , 



where the last passage is due to the change of variable 
t~ a = u. Let us now use, as in the previous appendix, 
the integral representation: 



-1—1/ a _ 



1 



r 



' a+l\ 



dzz 1 ' 



Through this transformation we arrive finally to the re- 
lation 



dt ( 1 — exp | / 
2 



t 



'a+l' 



dz ■ 



\t\ a + l 

z -l+l/a 

1 + z 2 



(35) 



Note that the last integral is diverging for a < 1/2, in- 
dicating that the problem is not well defined for these 
values of a as F is not a well defined stochastic quantity. 
This means that the sum in Eq. I|33|) needs a L depen- 
dent normalization to become a well defined stochastic 
variable. In fact, differently to the shuffled lattice case, 
where the typical mass fluctuation on regions of size R is 
proportional to R°, in the Poisson case this is due to the 
fact that such fluctuation is proportional to i? 1 / 2 . The 
field due to the mass fluctuation in a sphere of radius R 
on the origin of the sphere is of order R~ a for the shuf- 
fled lattice and R- a + 1 / 2 i n the Poisson point process. 
This explains why for a shuffled lattice the problem is 
well defined for any a > and not only for a > 1/2. 
This also says that for a < 1/2, in order to have a well 
defined stochastic field also in the Poisson case, we have 
to divide the field in Eq. |JSJ by L~ a+1 / 2 where L is 
the system size. The same argument can be used in d 
dimensions, for which the same mass fluctuations are re- 
spectively proportional to i?( d_1 )/ 2 in the shuffled lattice, 
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and R d / 2 in the Poisson case. This says that the prob- 
lem is well defined, without L dependent normalization 
of the field, for a > (d— l)/2 for the shuffled lattice and 
a > d/2 for the Poisson case (in the gravitational case in 
d = 3 faced by Chandrasekhar we have a — 2 and the 
field is a well defined stochastic quantity in both cases). 

Let us now go back to Eq. (|35[1 for a > 1/2. By using 
Eq. I|29|) we can rewrite it as 



dt ( 1 — exp i 



oF 



which is exactly of the form of the characteristic function 
of a Levy stable PDF [17| • By expanding Eq. (|3Tj|) to the 
first non-vanishing order larger than zero, we have: 



W P {k) 



«r(^)sin(£) 1 ; 



If we now invert this Fourier transform as explained in 
the preceding appendix, we can conclude that 

W P {F) ~ BF- {a+1)/a for large F , 



where we have also used cos[7r(l/a— l)/2] = sin[7r/(2a)]. 
Consequently 



Wp(k) = exp 



np 



«r(^)sin(^) 



(36) 



with 



B 



(37) 
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